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Abstract: In many control applications, the sensor technology used for the measurement of 
the variable to be controlled is not able to maintain a restricted sampling period. In this 
context, the assumption of regular and uniform sampling pattern is questionable. Moreover, 
if the control action updating can be faster than the output measurement frequency in order 
to fulfill the proposed closed loop behavior, the solution is usually a multirate controller. 
There are some known aspects to be careful of when a multirate system (MR) is going to 
be designed. The proper multiplicity between input-output sampling periods, the proper 
controller structure, the existence of ripples and others issues need to be considered. A 
useful way to save time and achieve good results is to have an assisted computer design 
tool. An interactive simulation tool to deal with MR seems to be the right solution. In this 
paper this kind of simulation application is presented. It allows an easy understanding of 
the performance degrading or improvement when changing the multirate sampling pattern 
parameters. The tool was developed using Sysquake, a Matlab-like language with fast 
execution and powerful graphic facilities. It can be delivered as an executable. In the paper 
a detailed explanation of MR treatment is also included and the design of four different 
MR controllers with flexible structure to be adapted to different schemes will also be 
presented. The Smith's predictor in these MR schemes is also explained, justified and used 
when time delays appear. Finally some interesting observations achieved using this 
interactive tool are included. 
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1. Introduction 

A multirate sampling (MR) system is defined as a hybrid system composed of continuous time 
elements, usually the plant, and some discrete time components, usually the controllers or the filters, 
where two or more variables are sampled or updated at different frequencies. It can be also considered 
that the discrete actions are not equally spaced on time and/or delayed. Moreover, in a great number of 
computer control applications the approximation of a regular pattern of sampled signals is assumed. 

A non-very restrictive assumption to simplify the treatment is to consider that the sampling pattern 
is periodic. That is, the process variables are sampled and/or updated at different and/or irregular 
intervals, but there is a global period T 0 with cyclic repetition. It may be also considered that there is a 
delay between the sampling and the updating of variables, but a global periodicity is still assumed. The 
case of asynchronous sampling/updating, with a random occurrence of the discrete actions, is much 
more complicated and it will not be considered in this paper. 

In a basic digital control system, a perfect uniform sampling and updating pattern of the involved 
variables is assumed, but it should be pointed out that, in practical applications, the synchronicity of 
the set of discrete actions is not perfect or it can be modified in order to improve the performances. 
Thus, MR is an important issue not only for research purposes but also from a practical point of view. 
MR may be present in a wide range of applications and the users must understand its consequences in 
an easy way. Chemical analyses, or samples obtained by artificial vision with post-processing 
requirements need a time interval that for a real-time process control request could be long. Some other 
similar problems appear when the sensors are spatially separated from the controller algorithm device: 
distillation columns, UAVs, network based control schemes, etc. 

The control target is to achieve similar performances to those the faster single rate controller would 
provide. However, in these cases, the theoretical analysis of the controlled system performances is 
much more computationally involved. The modeling, analysis and design steps can consume a great 
amount of engineering time. 

In order to analyze and study the different characteristics of the dynamic behavior of a MR, 
it is common to use time and frequency techniques and tools. These will provide a complete and 
global picture of the system behavior, showing up the interrelation among the different controlled 
plant performances. 

The combined use of control system design tools and dynamic system simulation, leads to the 
computer aided control system design environments, simplifying the design task. Additionally the 
possibility of simultaneous visualization in various windows of the effects in different performances of 
some design parameter changes helps to observe with more flexibility the change gradient over the 
system [1]. This facility provides the understanding of the usual steps in a design procedure. The 
perception of synthesis and analysis phases is simultaneous with the consequent effort saving 
with relation to classical simulation environments. In general, the complexity of the theoretical 
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developments justifies the use of interactive simulation techniques that also allow for acting over a 
high number of parameters with hard crossed relations. The global and simultaneous dynamic 
visualization of different kinds of time and frequency diagrams allows grasping a clear understanding 
about the effects of the concerned topic [2]. 

In this sense some years ago, Astrom and colleagues at the Lund Institute introduced some valuable 
concepts for control education task aid. In this context the significance of concepts like dynamic 
pictures and virtual interactivity must be highlighted. This original idea was implemented in packages 
as Ictools and CCSdemo, that Johansson et al. [3] and Wittenmark et al. [4], developed at the 
Department of Automatic Control at the Lund Institute of Technology, and Sysquake, developed at the 
Institut d'Automatique of the Federal Polytechnic School of Laussanne by Piguet [5,6]. The old use of 
computer aided control systems design was definitively improved. The dynamic picture allowed one to 
handle with the mouse a set of different nature graphic windows with some common parameter/s 
among them. Some change in a parameter manipulated by the user implied the fast — practically 
immediate — modified visualization in graphics influenced by that object. One of the main advantages 
is that the user does not need the implementation of code sentences. The complete effort is leading to 
testing and understanding of the system control ideas and principles that the application involves. 

In the MR case this kind of application appears indispensable. Some key concepts in order to 
model, analyze and design MR systems are overcome by the use of this interactive application. The 
working principles of MR are easily understood using this procedure. A specific Sysquake application 
was implemented for MR. This tool takes advantages of the fast execution and excellent graphic 
features that the use of Sysquake provides. In summary, the main motivations in writing of this paper 
are the following: (1) to study, in a very direct way, how the modification of the different sampling 
periods that are involved in a MR system can lead to unexpected behaviours; (2) to present in a unified 
way the performance of different MR controllers in order to choose by means of this tool the most 
appropriate design in each case; (3) to provide a tool that is not only useful to do research on this kind 
of systems, but also to improve the learning of MR concepts for beginners. 

The structure of the paper is as follows: the next section introduces some preliminary material and 
foundations, definitions and notation; Section 3 presents the problem statement and provides the 
detailed design procedure of the assumed controllers. Section 4 describes the developed application 
using interactive simulation techniques, and then some interesting examples derived with this tool are 
given in Section 5. A conclusion section closes the paper. 

2. Theoretical Foundations 

The scope and purpose of this work has been exposed in the previous section. In this section, the 
general problem, basic notations and operations among signals and processes are going to be 
introduced. After exposing the kind of problems that a practitioner finds when consider this topic, the 
elemental signal change of frequency operations and its properties are presented. Another subsection is 
devoted to the notations in process transfer functions in the MR topic, some elemental transformations 
between polynomials as well as the available relations between fast- skipped and slow or 
slow-expanded and fast signals. Finally the discrete lifting, traditionally introduced in an internal 
representation way, is adapted to our algebra. It is a section that is a survey to follow the design 
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procedure in Section 3. First of all, it must be noted that the systems this contribution deals with are 
known as MultiMate Systems, that is, systems where there are sampled or discrete signals referred to 
two or more different frequencies. An initial scheme could help to understand different issues related 
to this kind of systems (see Figure 1). 



MT 



U(s) 



Figure 1. An initial MR System. 
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One option in order to describe the different signals and systems in these environments is to use 
notation with superscripts. The signal (or system, when it is the case) Y T denotes either the Z-transform 
of the sequence yikT) derived from the sampling with period T of the continuous signal y(t): 



Y T ±Z{y(k)} = ^Ty(kT)z 



k=0 



(i) 



or the sampling rate transformation of a discrete signal Y (as will be explained below). With respect to 
Figure 1: 



Y NT =[G(s)U m \ =g nt \_u mt \ 



(2) 



where G NT represents the continuous process discretization (usually ZOH-discretization) at period NT: 



G NT =Z 



l-e 



-NTs 



G(s) 



(3) 



This single example enables one to understand that the sampling period transformation between 
discrete signals or the sampling operations involving blocks of different nature is quite common in MR 
systems. One way to achieve a proper handling of this kind of systems requires management of the 
great common divisor (gcd) and the least common multiple (1cm) of every sampling period occurring 
in the studied MR scheme. With these magnitudes, every sampling period is going to be repeated an 
integer multiple of times in one 1cm, 7b. There also will be a base period (gcd), T B such that 7 0 = PT B , 

being P an integer greater than one. With respect to Figure 1, if it is called To the 1cm period and 
B = 1cm (M,N), T B = To/B. 

Once the general MR problem has been introduced, the multirate input control case (MRIC), that is 
when the output sampling period is strictly greater than the input sampling period, will be assumed in 
the rest of this work. In these conditions, it is adequate to assume two sampling periods, input T 
and output NT, being TV a positive integer number greater than 7 if a MR system is studied (N is known 
as multiplicity). 



2.7. Signals: Basic MR Operations and Properties 



Some operators adapted for the transformation between periods with integer multiplicity must be 
introduced. In this sense, it is useful to consider basic operations such as "skip" and "expand", 
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introduced by [7], that look at the data in the same way as the classical downsampling and upsampling 
operations used in signal processing. 



If the Z- transform referred to period NT is defined as: 

00 

Y m ^Z NT {y{k)} = Y J y(kNT)z 



k=0 



(4) 



it is possible to formally express: 

• the expand (upsampling) operator creates a T-sequence from a ATT-sequence, as follows: 



[Y NT (z N )J=Y T (z N )=±y(kT)z 



k=0 



_ m [y{kT) = y{kT);\/k = m 
\ y(kT) = 0;Vk*AN 



the skip (downsampling) operator creates a A^T-sequence from a T-sequence, as follows: 
[Y T (z )f T =Y NT (z N )=±y(kNT)z- kN =Y NT (z N ) 



k=0 



For a graphical explanation of this kind of operators, see Figures 2 and 3. 

Figure 2. Expand operation. Case N=3. 



| 5 
> 

Q_ 

E 4 
cr 

3 



Output 
*■ Input 



(5a) 



(5b) 



0 1 2 3 4 5 



Figure 3. Skip operation. Case N = 3 
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The skip operation applied to the Z-transform of a signal can be obtained using the expression (6), 
due to Sklansky [8]: 



r t ~\NT 1 ^— \ i — 



z -Z N 



Some known skip-expand properties usually considered in this work are: 

[x nt (z n )Y nt (z n )J=[x NT J[y NT J(z n ) 
[x t (z)[y nt1 [(z)\ T =[x T Y T Y NT (z n ) 



(6) 

(7a) 
(7b) 

(7c) 



that is, (a) the skip does not commute, and (b) the expand commutes. The third is a clear rule with an 
easy proof that will be used. 

As an example guide to demonstrate every other property, the pattern "skip operation does not 
commute", (7a), is proved. Without loss of generality N = 2 has been considered: 



[x T Y T f=< 



(8) 



= \_{x 0 + x x z~ l +x 2 z~ 2 +x 3 z~ 3 + x 4 z~ 4 + ...}* {.y 0 + y x z~ l + y 2 z~ 2 + y 3 z~ 3 + j 4 z" 4 + ...}] = 

= [{^o) ; o + ( x Jo +^o) ; i)^ 1 + (x 0 y 2 + x Ji +x iy*)z 2 +(^o) ; 3 + x Ji + + x ,y*)z* + 

(^0^4 + ^3 +^2^2 +X ^l + X J^ + -f T = Wo + ( X 03 7 2 + X Jl + Wo)*'* + 

+{x 0 y A +x l y 3 +x 2 y 2 +x 3 y l +x A y Q )z~ A +... 

* [X r J" [F r ] 2r = {x 0 + Oz" 1 + x 2 z~ 2 + 0z~ 3 + x 4 z~ 4 +...}*{ y 0 + Oz" 1 + y 2 z~ 2 + 0z~ 3 + j 4 z" 4 +...} = 

= {^D +0 ^ 1 +(^2^0 + ^03 7 2 )^" 2 +0 ^ 3 + (^ 4 3 ; 0 +V4 +^ 2 J 2 )Z" 4 + ...} 

Once some basic operations have been introduced, it is also very interesting to expose the 
application over transfer functions. From here, the expressions will be managed without the subindex 
"1" for the variable z (r-period). Every formula can be assumed in z referred to T, or z N = z N with 

respect to NT; the user knows the argument of every signal. 

2.2. Process Transfer Functions: Notations and Transformations 

Assuming the continuous time process in Figure 1, for any pair of sequences like those in that 
Figure, which are respectively considered as process output and input, a discrete time (DT) transfer 
function of the process plus the hold device can be written: 

1. The fast sampling DT (FSDT) model is defined by: 



U T (z ) 



= G T (z ) = [H T G p ( S )f = 



(z) l +YL a * z~ { 



(9) 
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where B T (z),A T (z) are polynomials in z Following the notation in (1), these polynomials also 
represent finite sequences with ^-elements b iT ,a iT , respectively, for i = 1, n, being aoj = 1, and 

fe 0 ,r = 0. The expand and skip operators, (5), can also be applied to polynomials in z~ l . 

2. For the same process, a slow sampling DT (SSDT) model can be similarly defined by: 

y t NT / N, =G ( Z ) = l H NT G p( S )\ =" 



(10) 



ttNT/ N\ x y L ^Vi pv 'J ANT/ \ „ 

C/ (z ) A (z„) l + X M a, OT z - 

where B NT (z N ), A NT (z N ) are polynomials in z^ 1 , and mathematically in z _// . The proper treatment 

expressed when defining the skip and expand signal operations will be assumed for these polynomials. 

When the operation [•J 1 is applied to a system it must be considered the Z transform of the impulse 
response of that system. 

The FSDT transfer function poles are denoted by oc i T . That is: 

^(^ = 1X1^^) (11) 
If the SSDT poles are a iNT , such that: 

A NT (z N ) = lljz N ^ NT ) (12) 

Note that, dealing with the same continuous time system, a t m = a { T N ; V/ = 1, . . . , n . 
Two technical assumptions are made: 

27tkj/ 

a. If a is a pole of G T (z), then ae ^ N J=1...(N - 7),is not a pole of G T (z) ; 

b. All poles are different, i.e. , a iT * a. T \i* j . 

Assumption (a) is required to avoid aliasing, [9]. By assumption (b) the notation is simplified. 
The following useful relationship is easily derived: 

K(z j.m^.^gt.^.rttfM^^o (13) 

W =l (z -a iT ) A KZ ) A (z) (= i 
The FSDT model may be also expressed by: 

Y T (z) =cT B T (z) = B T (zWl(z) = B T (z) = YHIAt z~' 

U T (z) Z A T (z) A T (zWl(z) [A NT (z N )J 1 + z~ iN ( } 

where G T (z) = Z T [H T (s)G p (s)] , being Z r the Z transform referred to a period T and H T a T zero 
order hold. 

In what follows, some important results are proved: 
From (14), it yields (every variable expressed in z): 

B T U T =[A NT jY T (15) 

If a skip operation is applied to the T-time sequences, that is, doing a NT re-sampling, the result is: 
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[B T U T J T =[[A NT J Y T ] ;[B T U T Y =[[A m J] [Y T J T ; [eFU^ ^Jy^ (16) 

where is A , but expressed by means of the z-variable. The physical meaning is the 

consideration of a NT sampler at the process output, that is, a slow output: 

[B T U T Y =A NT Y NT (17) 

['~ T t ~\nt 
B U J , because the skip operation 

does not commute, (7a). From this fact is derived that it is not feasible to plan a transfer function 
between a skipped fast input and a slow output. 

The opposite situation is viable: the transformation of a slow frequency DT sequence into a fast 
frequency DT sequence. The dual rate zero order hold (DRZOH) device is defined by: 



u 



[u NT J 



1-e 



-NTS 



= = (1 + z* + + ) =: W- 1 (z) =: fc 1 J (18) 

1 - z 



where the output is defined by U T 0 . 

In this case, it is possible to obtain a transfer function of the process plus the DRZOH device: 

Y T =[jy OT G p (j)t/ OT ] r =[H NT G P (s)f [U^J 

because, as it is known, the expand operation commutes (7b). Also, using (18): 

[H^Gpis)] 7 =^f [H T G P (s)f 

Thus, a dual rate discrete time (DRDT) operator is defined by: 



(19) 



(20) 



yT,NT 



A T W T A 



(21) 



TNT 

G ' describes the transfer function from an expanded slow input (NT) to a fast output (T). 

Using a similar notation, the DRZOH operation, (18), can be represented by h t,nt . Clearly, it is 
also verified that: 



NT 





NT 


Y T 





B 



NT 



NT jjNT 



NT 



, NT 



(22) 



that is the SSDT model. 
2.3. Discrete Lifting 



When facing the modeling step of a MR system, most authors suggest the use of the so-called 
"Lifting" or "Discrete Lifting" method [10,11]. With this procedure every signal is referred to the 1cm 
of the periods of the MR system, and consequently it is "lifted" at 1cm period. This procedure can also 
be explained using the skip-expand operators. With the intention of showing the application of Lifting, 
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a general MR scheme where NT is the 1cm is assumed. One T-discrete signal Y T will be modeled in 
the lifted field by: 



y 0 + y N z +■■■ 
yi + y N+ iZ~ N +■■■ 



y N - l +y2 N - l z~ +■ 



(23) 



That is, every sequence y u (z N ) is obtained expanding N previously skipped elements from the 

T-original signal. 

In this way it is possible to say that: 

Y(z) = (l z x z~ 2 ... 

because: 



-(N- 



(24) 



y(z) = X y k z k = ( y 0 + y N z~ N +•••)+ z~ l ( x + y N+ ^~ N +...)+...+ z- {N - Y) (y N _ x + y 2N _ x z~ N +...) 



k=0 



(25) 



It must be noted that T is the gcd period of the scheme. 

Anyway, there are different options for Lifting application. All of them assume the technique of 
"Vector Switch Decomposition" [12]. That is, a multivariable system is achieved. It is usual to 
consider a state-space transformation, but one problem is derived of this. It is usual to find a mistake in 
some works that assume state-space control techniques; in the MR system there are not just variables 
with different dimensions, there are variables considered in different sampling times on one 1cm period 
(sometimes called meta period or frame-period). That is the reason why in this contribution only the 
external representation will be considered [8,13-15]. 

2.4. Modeling MR Scheme at Fast Period 

As exposed in Section 2.3, it is viable to express a fast signal as a sum of N slow signals. This result 
is especially interesting in order to solve the problem exposed in (17). In fact, using (7c): 

[G T H T f T =lG^ + z- l Gr+... + z- (N - l K T jH T Y = 
[G 0 NT J H T + z~ l [or J H T + . . . + z - (N ~ l) [g n n t _, ]H T \ T = 

G 0 NT [H T f T + Gr [z~ l H T f T +... + [z^H T f r 



(26) 



As it is proved, the terms can be separated. Actually, it is a laborious procedure but it assures the 
feasibility to get a fast sampling period modeling from a dual-rate closed loop, as it will be shown in 
Sections 3 and 4. 
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3. Problem Statement 

The problem this contribution deals with appears when a computer control application, where the 
control algorithm implementation should be a T single rate control, is not viable. The restricted 
frequency sensor is the main reason. In different fields like chemical industries, artificial vision, 
network-based and distributed installations, the sensors need a certain amount of time (for computation 
requirements or due to sensor geographic location) that makes the ideal sampling period unfeasible. 
This is the environment where a MR control can be a valuable option. In order to overcome this 
problem, the idea is to reach the same T-behavior but measuring the controlled variable N times 
slower. As it easy to understand, a non-conventional controller is needed; the controller uses a slow 
input but it must deliver a fast output. The basic scheme is showed in Figure 4, where a dual-rate 
control is introduced. The plant is represented by an n-order single-input-single-output LTI continuous 
system (CT), with transfer function G p (s). 

Figure 4. Dual rate controller structure. 
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Based in this MRIC structure, that is slow measurement and fast control updating, in [8] a new non- 
conventional structure composed by slow and fast parts for the dual rate controller was exposed; an 
expand operation is required to assure the composition of two different frequency elements. The rest of 
the closed loop performs as follows: the controller output is updated at a period T through the fast hold 
device, H T , and the system output, y(t), is measured at period NT, being represented by a fast sampler 
followed by an TV-sampler skip operation, and compared to the reference R(t) sampled at the slow rate, 



R . The dual rate controller, G MR , processes the error at slow rate, E , and generates N fast control 

actions, U T , inside the meta period NT. It must be noted that just the case where in all blocks input and 
output sampling periods are integer is considered; that is, the more complex case when rational ratio 
appears is out of our scope. 

As Figure 4 shows, and using the skip-expand properties (7), the controller output can be given by: 



u t =o T 2 [urJ =Gi[ G r][E NT ] =Gi[ G r][R NT -[/J'] 



therefore: 



and, since y 



T 

y dr 



DR 



= GIH t G t 2 [g^][r nt -[yl R f] 
\ y T DR f ' ^ e Equation (28) cannot be expressed in a closed form. 



(27) 



(28) 
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Nevertheless, if the dual rate modelling is considered at slow rate: 



(29) 



and finally: 

bog J 

R NT i + [g t p h t g t 2 ] nt g! 

With a similar procedure, the following expression is reached: 



NT 



(30) 



T 

y dr 



G T H T Gl[G^ 



[r nt ] iJ[g t p h t g t 2 J t ] [or] 



(31) 



In order to avoid the denominator complexity; (26) could be used. 

These expressions are useful when the simulation tool is implemented. Once the dual-rate closed 
loop modeling has been established (30), the design is faced. The goal is to force the skipped fast 
output, [l^J^, to be the same than the slow single rate control loop output, Y , that is, [flf r /f r J Vr 
should match to M NT R NT . One of the options is to consider a model based controller. In this way, if 
M(s) is considered as the reference model for the controlled system and M NT and M T are the ZOH 
discrete equivalent transfer functions from M{s), according to [8] the dual rate controller will be 
formed by: 



a fast part given by: 



G T 2 {z ) = 



M T (z) 



(32) 



a slow part given by: 



cr(z N )= 



l-M NT (z N ) 



(33) 



a rate converter with the form: 



H T (z ) = 



R T (z) 
[R NT J(z) 



(34) 



It is easy to observe that the dual rate controller is tuned for one type of command signal R(t) , 
since the rate converter (34) depends on the selected reference. In this paper this particular kind of rate 
converter has been used, but others can be used. 

3.1. Model Based Controller 

In this case, the design procedure is based on a continuous closed loop including a controller G R (s) 
(usually, a PID-type controller). From this controller a closed loop transfer function M(s) is obtained. 
Then, ZOH-discretization assuming periods T and NT must be computed in order to obtain M T and 
M NT , and respectively. So, applying (32)-(34) the design step is completed. 
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If the process is non-minimum phase, the cancellation of unstable pole-zero pairs must be avoided. 
Thus, the fast part of the controller could be alternatively computed by (35): 

Gl (Z\ ) = — ~ — — = R ^ — (35) 

G T ( Zl ) l + G T R (z)G T (z) 

If the slow part (33) is conserved, the design method is obviously not the introduced one, that is, the 
output does not match that predicted by the closed loop transfer function. 

3.2. Cancellation Controllers 

From the continuous time closed loop, it is possible to follow another design method. If some 
desired M(s) is considered (note that now M(s) is not derived from a closed loop containing a specific 
controller and plant), the expressions (32)-(34) give again the dual-rate controller. Obviously, it is a 
classical cancellation-type controller, and therefore the basic rules to select a proper M(s) must be 
assumed. For more information about different kinds of multi-rate cancellation controllers (Minimum 
Time, Finite Time), reader is referred to [8] and [16]. 

3.3. PID Dual Rate Controllers 

With the proposed structure (28)-(30), another alternative way to obtain the slow and fast controller 
parts can be starting from a PID controller but without reference model. When a PID-type controller is 
considered, it is logical to think with the derivative action working in high frequencies and the integral 
action in the lower ones. Therefore, a possible decomposition according to classical discretization 
methods, [17] could be defined by: 

I T i J (36) 



G pi(Zn) = K pi - f 



V T j 



(37) 



G PD (Z) = K PD 

z 

For more information, reader is referred to [18]. 
3.4. RST Controllers 

The RST controller [19] is a two-degree-of-freedom controller. Its design is an input-output 
model-based pole placement design procedure, requiring the resolution of a Diophantine equation. In 
this paper (as shown in Figure 5) to avoid confusion, the RST polynomials will be named as RST 
(since R is used for the reference signal, and T for the sampling period). The polynomials r,s are 
deduced when solving the Diophantine equation. The polynomial f is designed to tune the control 
system gain, and to avoid unstable zero cancellations. 
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Figure 5. Block diagram for the control system: multi-rate controller including an 
RST stage. 
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In this case, RST control law is designed (see [1]) at period NT, and yields: 



~NT T 



NT 



NT 



(z N ) 



R NT (z N ) R NT (z N ) 



[r T J T 



(38) 



This slow frequency control signal is modified by a dual rate part obtained from (32)-(34). 
So, from: 

yNT r^NT 



M NT = RST 



1 . r^NT s^NT 



(39) 



the slow side of the multi-rate controller can be rewritten in this way: 

G NT 



1 



l-M 



NT 



= G 



NT 



RST 



Thus, the multi-rate controller can be redefined as follows: 



-yT,NT 
7 MR 



M 



R 1 



That is using (13): 



M 



NT 



[ RNT \ 



T [b nt J 
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NT 
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NT 
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In addition, G\ m (z) can be represented by: 



G 2rst ~ 



M 



Leading to: 



2rst 



B T M A T [B NT J AI(zW t Ma B t m W t Ma [B NT ] 



K B T [ B "/J A t (z)W t a ~ B T W T A [ B »/] 7 



(40) 



(41) 



(42) 



(43) 



(44) 



which define the slow and fast parts of the dual rate controller. For more information, the reader is 
referred to [20]. 
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3.5. Processes with Time Delay 



In the MR control design field it is difficult to find contributions about processes with time-delay. 
The projected interactive simulation application tackles this problem, even in the general case when 
the delay is integer or non-integer with respect to the fast T and slow NT sampling periods in the 
dual-rate scheme. In general the delay, d, will be: 

d = LxNT + JxT + rn (45) 

being L and / the integer multiplicities with respect to the slow and fast sampling periods. 

Revisiting the control scheme already shown in Figure 4, in this case the Smith's Predictor (SP) 
proposal is adapted. Basically, this solution looks like a classical SP scheme, but now the novelty is 
provided by the algebra introduced in Section 2. In general, if d ^ 0 , the ZOH discretization of the 
process must be faced using the modified Z- transform [21]. 

With respect to the Figure 4, if the expression (28) is revisited and a skip operation is assumed: 



? NT 



E NT = R NT 

If the process includes a transport delay d, then: 



\y t T 



(46) 
(47) 



[Y T J T = (Z0H T G p (s) e ds J H T G T 2 
z^z' 1 (Z0H T G p (s) e- ns ) T H T G T 2 



s~iNT rpNT -L 

Lr { h — Z 



(Z0H T G p (s) e -(^^^ J h t G t 2 
z' J (Z0H T G p (s) e ms J H T G\ 



g nt e > 



qNT gNT 



(48) 



That is, the branch to subtract in order to avoid the effect of the delay. 

It must be noted that {zOH T G p {s) e~ ms ) represents the process with T-zero order hold Z-modified 

transform. In Figure 6, G T represents the T-ZOH discretization with a delay (/xl + m); obviously 

if m ^ 0 , a T-modified Z transform must be considered. This procedure can be directly applied to the 
cancellation, model-based and PID controllers. In the case of RST, a special figure must be adapted 
(see Figure 7), but the proof is very similar. 

Figure 6. Dual rate scheme with Smith's Predictor. 
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Figure 7. RST Controller in a loop with dead time. 
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4. Description of the Application 

To better understand the multi-rate control aspects previously presented, an interactive simulation 
application has been developed. An executable version can be downloaded online from [22]. The main 
goal of this section is to detail the main features of this application, which has been programmed using 
an interactive simulation tool named SysQuake [6]. SysQuake uses a Matlab-like programing 
language, which is provided with some specific commands in order to perform a high level of 
interactivity in the developed application. This fact enables students to more easily learn about 
advanced, complex control concepts and techniques, and researchers to exploit and improve their 
achievements [23-25]. 

Being aware interactivity is difficult to be shown in a writing text, next the different sections of the 
interactive application will be described. The application user interface presents a main window which 
is split into two parts (see in Figure 8): parameter section (upper left-hand part) and graphic section 
(lower left-hand part, and right-hand part). In addition, at the top side, the interface shows a menu bar 
and a toolbar, and at the bottom side, a status bar. These bars are provided by SysQuake and 
their working mode is standard, except for the Settings menu, which can be intentionally defined by 
the programmer. 



Settings menu: by means of the option Controller of this menu, one of the four multi-rate 
controllers presented in the previous sections (PID controller, Cancellation controller, RST 
controller, or Model-based controller) can be chosen. Moreover, increments/decrements for 
every slider included in the interface can be defined by means of the option Slider Config of 
this menu, which provides two variation modes for the sliders: by percentage or by fixed value. 
Parameter section: this section (located in the upper left-hand part; see Figure 8) enables to 
introduce the input data required in each case. The section is split into two parts in order to 
consider input data of different nature: process parameters (right-hand part) and controller 
parameters (left-hand part). 
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Figure 8. Application user interface for the PID controller. 
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Regarding the process parameters, two radio buttons enable selection of the process model 
complexity (first or second order). Data entries to define the consequent process parameters can be 
carried out via sliders or text boxes. Below these entries, a checkbox named Integrator is located, 
which can be activated to add an integrator to the selected model. In addition, the sampling time T 9 the 
multiplicity N, and the process time-delay D can be defined via slider (the first one) or text boxes (the 
last ones). While inserting every process parameter, the current process model is presented in blue text 
as a continuous-time transfer function G p (s) (below the checkbox). 

With respect to the controller parameters, the input data depend on the chosen controller. Thus, 
firstly, for PID controllers (illustrated in Figure 8), four radio buttons allow choosing the controller 
complexity (P, PI, PD, or PID), and consequently, the data entries for each basic control action can be 
carried out via sliders or text boxes. While inserting every parameter, the current continuous-time 
transfer function for the controller G r (s) appears in black text below the four radio buttons, whereas the 
closed-loop transfer function M(s) appears in red text above the radio buttons. Secondly, for 
Model-based controllers (see Figure 9), their interface is similar to the previous one, but, since only a 
PID controller structure can be treated, the four radio buttons are eliminated. Another difference is the 
existence of a checkbox named Oscillations, which can be deactivated to eliminate process output 
oscillations (remember (35)). Finally, for Cancellation and RST controllers (see Figure 10), a similar 
interface to that described when inserting the process parameters is used, but now with the goal of 
defining the desired continuous-time reference model M(s)(which, as previously, appears in red text as 
a transfer function). As a summary, Table 1 gathers every parameter employed in the developed 
application. 
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Figure 9. Application user interface for the Model-based (PID) controller. 



ysquake 5 Pro - [Multi-rate Controllers. sq] 
) File Edit Settings Figure Layout View Window Help 



-EBB 



Model Based Controller 
PID 



kp: 
Ti: 
Td: 



w ■ (s + 0.31) (52 + 4.09s + 12 25) 

Oscilations: J 



-6 

Switch to Root L< 



Process Output 



o Manual ! AutoAdjust 



1st Order * 2nd Order 



Wn: 0.8 
Integrator 




Figure 10. Application user interface for the Cancellation controller. 
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Table 1. Parameters available in the interactive simulation application. 


Parameter 


Description 


kp, Ti, Td 


PID proportional gain, integral time, and derivative time, respectively. 


P, PI, PD, PID 


Four radio buttons to select the PID-type controller complexity. 


K, T.C. 


Static gain and time constant, respectively, for a first order model/process. 


K, D.F., Wn 


Static gain, damping factor, and natural frequency, respectively, for a second 


order model/process. 
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Table I. Cont. 



Parameter 



Description 



1st Order, 2nd Order 
Integrator 



Two radio buttons to select the model/process complexity. 
A checkbox to add an integrator to the process (when activated). 
A checkbox to eliminate process output oscillations in Model-based controllers 



Oscillations 



(when deactivated). 
Multiplicity. 
Sampling time. 
Process time delay. 
A link to show the discrete- time root locus plots. 
A link to come back to the continuous-time poles and zeroes map. 
Two radio buttons to select operation mode for figure scales. 



N 



T 



D 



Switch to Root Locus 
Switch to Pole-Zero Map 
Manual, AutoAdjust 



Predictor mrate 



Predictor slow 



Predictor fast 



A checkbox to include a Smith's Predictor in the slow-rate control algorithm 

(when activated). 

A checkbox to include a Smith's Predictor in the fast-rate control algorithm 

(when activated). 

A checkbox to include a Smith's Predictor in the multi-rate control algorithm 

(when activated). 



• Graphic section: controlled system output (named Process Output), control action signal 
(named Controller Output), and poles and zeroes map (named Poles & Zeros) are the figures 
plotted in this section in the main window. In addition, for PID and Model-based controllers, an 
option named Switch to Root Locus is available to replace the Poles & Zeros figure with the 
slow-rate and dual-rate root locus (named Root Locus Mono and Root Locus Bif, respectively; 
see Figure 11). When switching to this plot, a new option named Switch to 
Pole-Zero Map appears and enables to come back to the previous figure. 

Regarding the Process Output and Controller Output figures, they show the simulation responses 
for the single-rate (slow-rate in green, and fast-rate in cyan) and the dual-rate case (in black), which are 
achieved by the chosen control algorithm when following a step reference. In these graphics, an 
impulse-like load disturbance can also be considered by dragging on the cyan circles to interactively 
define its amplitude and application time (at instant multiple of NT). Moreover, scales can be redefined 
by dragging or clicking on the black diamonds. Another graphic feature can be observed in Figure 11: 
if the mouse is positioned on the process output or on the controller output, a text appears indicating 
the time instant and the output value for this mouse position. 

With respect to Poles & Zeros figure, it shows the location of the different continuous-time poles 
(crosses) and zeroes (circles) for the process (in blue), for the controller (in black; if PID or 
Model-based controllers are chosen), and for the closed-loop model (in red). Note that the colors are 
consequent with those used for the different transfer functions. By dragging on poles and zeroes the 
related continuous-time transfer functions and parameters change interactively, and vice versa. Finally, 
Root Locus figures show the discrete-time root locus for the slow-rate and dual-rate cases (observe that 
in this plot closed-loop poles appear and are represented by means of black triangles). Contrarily to the 
previous figure, poles and zeroes cannot be now dragged, but they change their location interactively 
when the appropriate parameters are modified. 
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Figure 11. Slow-rate and dual-rate root locus graphic. 
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5. Examples 

In this section an example for each controller will be presented. Output responses for the single -rate 
and multi-rate cases will be compared and analyzed, focusing on multi-rate control benefits. 

5.7. Model-Based Controller 

Let us consider the example shown in Figure 9, where the process is defined by: 



G D (s) = - 



1.5 



-0.5s 



(49) 



0 + 0.5)0 + 1.5) 

An acceptable continuous-time PID controller is given by K =8, T D =0.2, T t =3.2 (considering 

i V 



u(t) = K P 



i r 

e(t) + T D e(t) + — ^ e{r)d 2 



), which yields: 

r ( . Q 0 + 5)0 + 0.31) 
G r (s) = 8 



(50) 



From (49)-(50) a continuous-time closed-loop transfer function M(s) can be calculated. Then, if the 
sampling time is T = 0.1 s, and the multiplicity is Af = 4, the dual-rate controller is designed from 
(32)-(34) (note that the single-rate controllers can be designed using (32) and (33) with the specific 
rate for each controller). Responses for each case are illustrated in Figure 9, where stable outputs are 
always obsered due to the use of the Predictor options to face the process time delay (D = 0.5). In this 
example, the dual-rate controller is able to achieve a response which is quite similar to that obtained by 
the fast-rate controller when following a step reference. However, when applying a load disturbance, 
the dual-rate output worsens with respect to the fast-rate one, since the dual-rate controller is designed 
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for step references and not for impulse inputs like the load disturbance is. Regarding the slow-rate 
output, it is clearly worse for the both cases (when following a step reference and when applying a load 
disturbance), since it shows around 10% higher overshoot and 2 times longer settling time with respect 
to the dual-rate output. In any case, every output shows oscillations (ripple) due to cancelling the 
process dynamics. Deactivating the checkbox Oscillations, the controllers are designed following (35), 
and then the ripple effect totally disappears (see Figure 12). Moreover, in this case, whereas the 
slow-rate output is on the verge of instability, the dual -rate output presents around 5% less overshoot 
than the fast-rate output when following a step reference. 

Figure 12. Application user interface for the Model-based (PID) controller with 
no oscillations. 
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5.2. Cancellation Controller 



Let us consider the example shown in Figure 10, where the process is defined by: 

0.47 



GAs) = — 
p (s + 2.24)0 + 0.16) 



-6s 



(51) 



In order to decrease around 10 times the settling time, the next closed-loop transfer function 
is chosen: 



M(s) = 



4.97 



(s + 3.39)O + 0.98) 



(52) 



If the sampling time is T = 2 s and the multiplicity is N = 4, as in the previous example a different 
controller can be designed using (32)-(34). Responses for each case are illustrated in Figure 10, where 
as a consequence of including the Predictor options in every control algorithm, every response is 
stable. In this example, the dual-rate controller is also able to achieve a better response (no overshoot 
and around 20% shorter settling time) than that obtained by the fast-rate controller when following a 
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step reference. However, as previously, when applying a load disturbance, the dual-rate output worsens 
with respect to the fast-rate one (around 10% higher overshoot, and around 15% longer settling time). 
Regarding the slow-rate output, it is clearly worse considering the settling time index (around 50% 
longer with respect to the dual-rate output). The overshoot index is very similar. 

5.3. PID Controller 

Let us consider the example shown in Figure 8, where the process is defined by: 

136.17 



GJs) = — 
p s(s + 24.55)0 + 3.26) 



(53) 



An acceptable continuous-time PID controller is given by K =0.58,7^ =0.3.59, 7} =0.52 



(considering u(t) = K P 



e(t) + T D e(t) + y^e{r)d 



). From (31)— (32), and assuming K v \ = 1 and = K v , 



and T = 0.02 s and TV = 4, the dual-rate controller is designed (note that the single-rate controllers can 
be designed using (36) and (37) with the specific rate for each controller). Responses for each case are 
illustrated in Figure 8. 

In this example, the dual-rate response reaches time indexes whose values are between those 
obtained by the single-rate controllers. If some process time delay were considered, for example 
D = 0.06, the responses would significantly worsen (concretely, the slow-rate one would become 
unstable), as depicted in Figure 13 (where every Predictor option is now deactivated). In addition, and 
as well-known, the different root locus could not be generated. Nevertheless, if the Predictor options 
were activated, every response would recover the previous behavior (see Figure 14). 

Figure 13. Application user interface for the PID controller with process time delay and 
no predictor. 
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Figure 14. Application user interface for the PID controller with process time delay 
and predictors. 
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5.4. RST Controller 

Let us consider the next example, where the process is defined by: 

128.62 



G P (S)- 



0 + 2.42) z +11.08 



2 4-1 1 f>K 2 (54) 



In order to eliminate the overshoot and slightly decrease the settling time, the next closed-loop 
transfer function is chosen: 

= 60.72 

S 0 + 26.12)0 + 3.02) ( 55 ) 

If the sampling time is T = 0.25 s and the multiplicity is TV = 2, from (42)-(43) the dual-rate 
controller is designed (the single-rate controllers can be designed from (41) considering the appropriate 
rate in each case). The different outputs are illustrated in Figure 15. Once again, the dual-rate 
controller is able to achieve a better response than the fast-rate one when following a step reference. 
Concretely, the settling time achieved is around twice as short as that obtained by the fast-rate 
controller. However, as previously, when applying a load disturbance, the dual-rate output worsens 
with respect to the fast-rate one, yielding around twice as much overshoot and settling time. Regarding 
the slow-rate output, it presents an overshoot clearly worse than the dual-rate one. The settling time is 
also worse when following a step reference but quite similar when applying a load disturbance. 
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Figure 15. Application user interface for the RST controller. 




6. Conclusions 

One solution for the problem of control schemes with slow sensors is to assume a MR system, 
considering a restricted slow measurement sampling but also a faster control updating, the use of an 
interactive simulation tool in order to study multirate systems is a feasible option to make proper 
decisions about the correct control design. Different non-conventional controller structures have been 
provided. From application some unexpected results are achieved and explained. Even for an expert in 
this field the tool appears absolutely essential and time saving. For a beginner student or researcher it 
is entirely necessary if the study of this kind of systems is needed/desired. 

Acknowledgments 

This work was supported in part by the Spanish Ministry of Economy and Competitiveness under 
Project DPI2012-31303. The authors J. Salt, A. Cuenca, are grateful to Grant TEC20 12-3 1506, from 
the Spanish Ministry of Education. The work of A. Cuenca was supported in part by the Spanish 
Ministerio de Economia under Grant DPI201 1-28507-C02-01. 

Conflicts of Interest 

The authors declare no conflict of interest. 



References 

1. Dormido, S. Control learning: Present and future. Annu. Rev. Control 2004, 28, 115-136. 

2. Guzman, J.L.; Dormido, S.; Berenguel, M. Interactivity in education: An experience in the 
automatic control field. Comput. Appl. Eng. Educ. 2013, 27, 360-371. 



Sensors 2014, 14 



4109 



o 

3. Johansson, M.; Gafvert, M.; Astrom, K.J. Interactive tools for education in automatic control. 
IEEE Control Syst. Mag. 1998, 18, 33-40. 

4. Wittenmark, B.; Haglund, H.; Johansson, M. Dynamic pictures and interactive learning. IEEE 
Control Syst. Mag. 1998, 18, 26-32. 

5. Piguet, Y.; Holmberg, U.; Longchamp, R. Instantaneous Performance Visualization for Graphical 
Control Design Methods. In Proceedings of the 14th IF AC World Congress, Beijing, China, 5-9 
July 1999. 

6. Piguet, Y. Sysquake 5: User Manual, Calerga: Lausanne, Switzerland, 2013. 

7. Coffey, T.C.; Williams, I.J. Stability analysis of multiloop, multirate sampled systems. AIAA J. 
Guid. Control Dyn. 1966, 4, 2178-2190. 

8. Salt, J.; Albertos, P. Model-based multirate controllers design. IEEE Trans. Control Sys. Technol. 
2005, 73, 988-997. 

9. In, H.; Zhang, C. A multirate digital controller for model matching. Automatica 1994, 30, 
1043-1050. 

10. Khargonekar, P.P.; Poola, K.; Tannenbaum, A. Robust control of linear time-invariant plants 
using periodic compensation. IEEE TAC 1985, 30, 1088-1096. 

11. Bamieh, B.; Pearson, J.B. A lifting technique for linear periodic systems with applications to 
sampled-data control. Syst. Control Lett. 1991, 77, 79-88. 

12. Kranc, G. Input-output analysis of multirate feedback systems. IRE Trans. Autom. Control 1957, 
3, 21-28. 

13. Cimino, M.; Pagilla, P.R. Design of linear time-invariant controllers for multirate systems. 
Automatica 2010, 46, 1315-1319. 

14. Cimino, M.; Pagilla, P.R. Conditions for the ripple-free response of multirate systems using linear 
time-invariant controllers. Syst. Control Lett. 2010, 59, 510-516. 

15. Truman, A.W. An analysis of sub-rate control systems. Int. J. Control 1996, 63, 365-392. 

16. Cuenca, A.; Salt, J.; Albertos, P. Implementation of algebraic controllers for non-conventional 
sampled-data systems. Real-Time Syst. 2007, 35, 59-89. 

17. Isermann, R. Digital Control Systems', Springer Verlag, Berlin, Germany, 1981. 

18. Salt, J.; Cuenca, A.; Casanova, V.; Mascaros, V. A PID Dual Rate Controller Implementation 
over a Networked Control System. In Proceedings of the IEEE CCA/CACSD/ISIC, Munich, 
Germany, 4-6 October 2006; pp. 1343-1349. 

o 

19. Astrom, K.; Wittenmark, B. Computer Controlled Systems: Theory and Design, 3rd ed.; Prentice 
Hall: Englewood Cliffs, NJ, USA, 1997. 

20. Cuenca, A.; Salt, J. RST controller design for a non-uniform multi-rate control system. /. Process 
Control 2012, 22, 1865-1877. 

21. Jury, E.I. Theory and Application of the Z-Transform Method', Wiley: New York, NY, USA, 
1964. 

22. A Multi-Rate Control interactive simulation application. Available online: http://personales.upv.es/ 
jsalt/MR.exe (accessed on 18th February 2014). 

23. Dormido, S. Control Learning: Present and Future. In Proceedings of the 15th Triennial World 
Congress of IFAC, Barcelona, Spain, 21-26 July 2002. 



Sensors 2014, 14 



4110 



24. Dormido, S.; Dormido-Canto, S.; Dormido-Canto, R.; Sanchez, J.; Duro, N. The role of 
interactivity in control learning. Int. J. Eng. Educ. 2005, 27, 1122-1133. 

25. Salt, J.; Albertos, P.; Dormido, S.; Cuenca, A. An Interactive Simulation Tool for the Study of 
Multirate Sampled Data Systems. In Proceedings of the Advances in Control Education, Oulu, 
Finland, 16-18 June 2003. 

© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.Org/licenses/by/3.0/). 



